E = 0.65
G = 0.2
Q1 = 100
Q2 = 10e6
d0 = 10
s0 = 10
T0 = 400
rs = 100
import numpy as np
def Update(width, height ):
    '''
    规定数字信息素按照先传播后蒸发的原则，当前栅格的吸引信息素更新公式如下

    '''
    lame = np.ones(width, height)
    k = np.ones(width, height)
    p = np.ones(width, height) * s0
    d = np.ones(width, height) * d0
    g = np.ones(width, height)
    for i in range(width):
        for j in range(height):
            p[i, j] = (1-E)*((1-G) * p[i, j] + lame[i, j] * k[i, j] * d[i, j] + g[i, j])
            sum = G * (p[])